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MicroRNA profiling reveals that miR-21, miR486 and 
miR-214are upregulated and involved in cell survival in 
Sezary syndrome 

MG Narducci\ D Arcelli\ MC Picchio\ C Lazzeri\ E Pagani\ F Sampogna\ E Scala\ P Fadda\ C Cristofoletti\ A Facchiano\ 
M Frontani\ A Monopoli\ M Ferracin^, M Negrini^, GA Lombardo\ E Caprini^ and G Russo*'^ 

Sezary syndrome (SS) is an incurable leukemic variant of cutaneous T-cell lymphoma and its pathogenesis is still unknown. 
Diagnosis/prognosis may strongly ameliorate the management of SS individuals. Here, we profiled the expression of 470 
microRNAs (mlRNAs) in a cohort of 22 SS patients, and we identified 45 mlRNAs differentially expressed between SS and 
controls. Using predictive analysis, a list of 19 mlRNAs, including miR-21, miR-21 4, miR-486, miR-18a, miR-342, miR-31 and let-7 
members were also found. Moreover, we defined a signature of 14 mlRNAs including again miR-21, potentially able to 
discriminate patients with unfavorable and favorable outcome. We validated our data for miR-21, miR-21 4 and miR-486 by qRT- 
PCR, including an additional set of array-independent SS cases. In addition, we also provide an in vitro evidence for a 
contribution of miR-21 4, miR-486 and miR-21 to apoptotic resistance of CTCL cell line. 
Cell Death and Disease (2011) 2, el 51; doi:10.1038/cddis.2011.32; published online 28 April 2011 
Subject Category: Cancer 



Sezary syndrome (SS) is a rare and aggressive leukemic 
variant of cutaneous T-cell lymphoma (CTCL) characterized 
by the presence of neoplastic lymphocytes named Sezary 
cells, in the skin, lymph nodes and peripheral bloodJ 
Consistently with the complex phenotype observed, 
SS is associated with a poor clinical outcome with an 
estimated 5-years survival of only 25%.^ Until now, no specific 
therapy is available for this disease.^ Although numerous 
efforts have been made to identify pathogenically relevant 
genes, the cause of SS is still unknown. With this aim, high 
throughput analyses have been recently done on whole 
genome and transcriptome of CTCL identifying regions of 
recurrent allelic imbalance and altered profiles of gene 
expression. "^"^ 

miRNAs control gene expression at the posttranscriptional 
level by causing mRNA degradation and/or repressing mRNA 
translation. They have a wide recognized role in development, 
differentiation, growth, apoptosis and stress responses, as 
well as in T-cell homeostasis.^"'' ^ A growing number of studies 
indicate that deregulated miRNAs have a key role in cancer 
development, acting both as oncogenes and as tumor 
suppressor genes. This evidence is also supported by the 
fact that more than 50% of miRNA-encoding loci reside in 
chromosomal regions unbalanced during tumorigenesis. 
Furthermore, miRNA expression profiles are able to predict 
disease status and clinical outcome, as well as tumor 
progression and therapeutic responses. ""^ 

Expression patterns of miRNAs and their role in the 
pathogenesis of SS have only recently been addressed. ""^ 



Here, we present a miRNA profile of 22 SS patients 
identifying a signature that distinguishes malignant from 
healthy lymphocytes, as well as a series of miRNAs 
associated with disease and clinical outcome, and an in vitro 
evidence for a contribution of miR-21 4, miR-486 and miR-21 
to apoptotic resistance of CTCL cell line. 



Results 

Differential miRNA expression between SS and normal T 
cells. To identify miRNA signature in SS, we profiled the 
expression of 28 SS samples (22 patients and 6 follow-up) 
and six healthy controls (HCs) by the means of microarray 
analysis containing 470 human miRNAs. To determine 
whether global miRNA profiling could distinguish molecular 
groups, we performed an unsupervised analysis using heat 
map package of bioconductor. As showed by dendrogram 
(Figure la and Supplementary Figure S1) miRNAs 
separated the HCs from SS samples. In particular, HCs 
(n = 6) were clearly grouped and statistically validated in a 
distinct sub-branch, whereas SS patients were sub-divided in 
four distinct clusters (correlation = 0.95). The first one, very 
close to normal controls, contained 4 SS samples and the 
other three clusters comprised 9, 6 and 9 SS samples, 
respectively. Four of six follow-up (mSS43_1, 30_1, 39_1 
and 33_1) were clustered together with the respective 
original specimen, whereas mSS38_1 and 38_2 clustered 
in a different group, indicating that a change might have 
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Figure 1 Analysis of miRNA expression profiles and survival curves in the SS cohort hierarchical clustering of the samples, showing the most significant mlRNAs 
(P<0.01) (patients in columns, mlRNAs in rows), (a) Dendrogram obtained using Pearson centered correlation distance metric and average linkage cluster method (for the 
correlated heat map see Supplementary Figure S1). (b) Heat map of the 45 mlRNAs differentially expressed between SS and HC. (c) Heat map of predictive analysis between 
HC and SS mlRNAs selected by class prediction analysis, (d) KM analysis performed on four clusters of SS patients identified by unsupervised hierarchical cluster analysis 
showing the corresponding MS time (Wilcoxon test P= 0.047). (e) KM analysis performed with cluster n.2 (green) versus cluster n. 1 , 3 and 4 showing an MS time of 38.84 and 
74.57 months, respectively (Wilcoxon test P= 0.011) 



occurred in these patients. Statistical analysis performed to 
validate the four subgroups of SS patients is shown in 
Supplementary Table S2. 

A class comparison analysis between HCs and SS lead us 
to identify the 45 most deregulated miRNAs (P< 0.01) able to 
discriminate the two groups; 21 miRNAs were upregulated 
and 24 downregulated. They are displayed in the heat map of 
Figure lb and listed in Table 1. 

By using a prediction analysis of BRB array tools (see 
Materials and Methods), we identified an additional list of 19 
miRNAs able to predict in which class (SS or HC) one sample 
could be enclosed relatively to its expression profile (Figure 1c 
and Supplementary Table S3). This list included miRNAs 
already highlighted among the 45 most deregulated miRNAs. 
Most of them were reported to function either as oncogenes 
like miR-21,^'^ miR-214,^^ miR-486^^ and miR-18a belonging 
to miR-17-92 cluster, or tumor suppressors like let-7 members 
and miR-SI."*^ 



miRNA expression profiling delineates also survival 
classes in SS. We next sought to determine whether the 
four groups of SS samples identified by unsupervised 
analysis were associated with prognosis. As clinical data, 
including those regarding the survival, were available for the 
SS patients enrolled in this study (see Supplementary Table 
SI), we used these four SS clusters to evaluate their survival 
by Kaplan-Meier (KM) estimator. As shown in Figure Id, we 
found that four groups of SS patients displayed four different 
survival curves, resulting in a statistical correlation (Wilcoxon 
test P= 0.047) with a median for survival (MS) time of 71, 
38.8, 77 and 57.7 months, respectively (compare matching 
colors in Figures la and d). These data indicate that cluster 
n.1 (yellow), the closest one to normal controls, and cluster 
n.3 (blue) represented the patients with longer survival (MS 
of 77 and 71, respectively), followed by cluster n.4 (violet) 
with a MS of 57.73. Conversely, cluster n.2 (green) showed a 
marked survival decrease with an MS of 38.84 (Figure Id). 
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Table 1 The 45 most significant miRNAs (P<0.01) identified by class 
comparison analysis performed between HCs an SS 
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Abbreviations: HCs, healthy controls; mlRNA, microRNA; SS, Sezary 
syndrome. 

In bold are indicated chromosomal regions involved in gains (+) and losses (-) 
in SS. 



To increase class numbers and simplify data analysis, we 
compared the survival time of SS patients belonging just to 
cluster n.2 (green) with all other patients. The resulting KM 
analysis delineated the cluster n.2 as the group of SS 
patients with the worst outcome, showing an MS of 38.84 
versus a group with MS of 74.57 (P= 0.011) (Figure 1e). 

We also asked which miRNAs could discriminate these two 
different survival groups. For this reason, we performed a class 
comparison between the patients belonging to cluster n.2 (mSS 
n.43, 30, 22, 32, 39 and 45) versus aW other patients (clusters 1 , 
3 and 4). In this way, we identified a list of 47 miRNAs, all 
upregulated in individuals with the worse outcome (Table 2). 

As we observed a significant association between miRNA 
expression profiles and SS survival, we tried to establish if a 



Table 2 miRNAs differentially expressed between SS patients with unfavorable 
and favorable outcome 



Unique id 


Chromosomal 


Fold- 


FDR 


Parametric 


location 


change 




p-value 


hsa-nniR-199a* 


19p13.2 


18.39 


6.42E-03 


8.08E-04 


hsa-nniR-214 


1q24.3 


13.70 


1.21E-02 


3.56E-03 


hsa-nniR-199a 


1q24.3 


11.33 


3.78E-03 


1 .40E-04 


hsa-nniR-146a 


5q33.3 


4.30 


5.60E-03 


3.07E-04 


hsa-nniR-29b 


7q32.3 


4.01 


1 .08E-02 


1 .90E-03 


hsa-nniR-142-3p 


17q22+ 


3.66 


7.96E-03 


1.17E-03 


hsa-nniR-486 


8p11.21 + 


3.64 


1 .06E-02 


1 .73E-03 


hsa-nniR-155 


21q21.2 


3.30 
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3.12 
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7.03E-03 


9.38E-04 
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2.13 


1.11E-02 


2.14E-03 


hsa-nniR-17-3p 


13q31.3 
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1 .24E-02 
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13q31.3 


2.00 


1 .43E-02 


4.96E-03 
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1.97 
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11q23.1 


1.89 


1 .31 E-02 


4.36E-03 
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17q22+ 


1.86 


1 .20E-02 


2.89E-03 


hsa-nniR-194 


1q41/11q13.1 


1.83 


1 .20E-02 


2.65E-03 


hsa-nniR-195 


17p13.1- 


1.82 


5.81 E-03 


6.89E-04 


hsa-mlR-ISIc 


19p13.13 


1.81 


1 .31 E-02 


4.31 E-03 


hsa-mlR-ISIc 


19p13.13 


1.80 


5.60E-03 


4.57E-04 


hsa-nniR-28 


3q28 


1.76 


1.21 E-02 


3.25E-03 


hsa-nniR-107 


10q22.31- 


1.74 


5.81 E-03 


6.51 E-04 


hsa-nniR-185 


22q11.21 


1.71 


1 .08E-02 


2,00E-903 


hsa-mlR-SOb 


8q24.22+ 


1.71 


1.21 E-02 


3.04E-03 


hsa-nniR-25 


7q22.1 


1.70 


5.60E-03 


3.42E-04 


hsa-nniR-331 


12q22 


1.69 


1.21 E-02 


3.30E-03 


hsa-let-7d 


9q22.32- 


1.61 


1 .08E-02 


1 .96E-03 
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7.96E-03 


1.18E-03 
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8q24.22+ 


1.57 


1 .20E-02 


2.75E-03 


hsa-nniR-191 


3p21.31 


1.56 


1.21 E-02 


3.40E-03 


hsa-nniR-148b 


12q13.13 


1.53 


1.21 E-02 


3.60E-03 


hsa-nniR-769-5p 


19q13.32 


1.40 


1 .96E-04 


2.90E-06 



Abbreviations: miRNA, microRNA; SS, Sezary syndrome. 

In bold are indicated chromosomal regions involved in gains ( + ) and losses (-) 

in SS. 



limited number of miRNAs might be used to build a predictive 
model for SS clinical outcome. To identify the most predictive 
miRNAs, we used the survival risk prediction algorithm 
implemented in BRB-ArrayTools^^ on SS samples described 
above. The methodological principles of this algorithm have 
been described.''^ In brief, high- and a low-risk survival groups 
were defined by a multivariate model based on gene 
expression levels contained in each gene signature, the Cox 
regression coefficient for each gene (supervised principal 
component method) and two covariates as sex and age.This 
multivariate model was used in a leave-one-out cross 
validation process to assign risk-group membership for 
clinical samples. Statistical significance of the survival groups 
was assessed by the log-rank test. The top-ranking miRNAs 
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Table 3 Survival risk prediction analysis based on mlRNA expression (a) and 
age and sex covariates (b) 



(a) 



miRNA ID 


Chromosome 


% CV 


P-value 




location 


Support 




hsa-miR-30d 


8q24.2 


100 


0.0004324 


hsa-miR-486 


8p11.21 


100 


0.0007376 


hsa-let-7d 


9q22.32 


100 


0.0010049 
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17q22 


100 


0.0018857 


hsa-nniR-93 
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100 
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hsa-nniR-422b 


5q33.1 


95.24 


0.0213468 


hsa-nniR-15a 


13q14.2 


90.48 


0.0221395 


hsa-miR-142-5p 


17q22 


90.48 


0.0251359 


hsa-nniR-18a 


13q31.3 


85.71 


0.0263155 


hsa-let-7a 


9q22.32 


80.95 


0.0272317 


hsa-miR-301 


17q22 


90.48 


0.0274629 


hsa-nniR-34b 


11q23 


90.48 


0.0281414 


hsa-let-7i 


12q14.1 


71.43 


0.029799 


hsa-nniR-181c 


19p13.13 


57.14 


0.0392069 


hsa-nniR-33 


22q13.2 


42.86 


0.0410688 


hsa-let-7f 


9q22.31 


57.14 


0.0433045 


hsa-nniR-17-3p 


13q31.3 


47.62 


0.0445715 



(b) 



Sample 


Survival time 


Censoring 


Predicted 




(months) 


indicator 


risk 






(0 = alive, 








1 =dead) 




mSS-04 


54 


1 


High 


mSS-40 


81 


0 


Low 


mSS-26 


71 




High 


mSS-49 


74 




Low 


mSS-30 


88 




Low 


mSS-22 


30 




High 


mSS-32 


40 




High 


mSS-39 


26 




High 


mSS-43 


38 




High 


mSS-45 


42 




High 


mSS-33 


66 




Low 


mSS-51 


68 




High 


mSS-27 


174 




Low 


mSS-02 


75 




Low 


mSS-23 


79 




Low 


mSS-01 


90 




Low 


mSS-25 


26 




High 


mSS-21 


44 




High 


mSS-38 


101 


0 


Low 


nnSS-36 


72 


0 


Low 


mSS-48 


48 


1 


High 



Abbreviation: miRNA, microRNA. 

^List of 23 genes selected by fitting Cox proportional hazards models to be the 
best risk classifiers (alpha = 0.05). 

In bold are indicated mlRNAs mapping on chromosomal regionsof gain or loss. 



selected by this analysis were 23 (P< 0.042). Among them, 
miR-30d and miR-486 appeared the most significantly risk- 
associated (P< 0.0007) (Table 3a). In good agreement with 
our previous KM survival analyses (Figures Id and e), we 
observed that five out six SS patients (83%) showing poor 
outcome were classified as high-risk patients using this 
prediction model. Moreover, we obtained a consistent 
classification into low-risk group for 9/15 SS patients (60%) 



previously associated with a favorable outcome (Figures la 
and d). The misclassification of 4/7 (mSS 48, 21, 25 and 51) 
remaining patients might be explained by the moderate risk of 
disease observed for these individuals (Figures la and d), as 
also indicated by their intermediate survival time showed in 
Table 3b. These data obviously require further validations in 
independent sets of samples in order to develop an accurate 
and unbiased classification profile that might be used to 
predict at which risk-class a future patient will be associated. 

With the aim to further restrict the list of prognostic miRNAs, 
we also intersected the list of the 47 miRNAs showed in 
Table 2 and the list of 23 miRNAs obtained by Risk prediction 
analysis (Table 3a). This approach highlighted 14 common 
miRNAs, such as let-7d, let-7f, let-7i, miR-107, miR-142-5p, 
miR-155, miR-16, miR17-3p, miR-181c, miR-18a, miR-21, 
miR-301, miR-30d and miR-486 (Figure 2). 

Validation of tlie miRNA signatures by qRT-PCR. To 

confirm the miRNA array data, we performed a qRT-PCR 
using CD4+ purified from PBMCs. We investigated the 
expression of miR-21, miR-21 4 and miR-486 on a total of 23 
samples, of which 1 1 were new cases (samples SS41 and 
from SS52 to SS 64). Of these, 10 are alive (24-36 months 
follow-up), whereas SS63 has 68 months of survival. 
(Figures 3a-c). We also validated miR-18a, miR-342, Iet7b 
and Iet7c in three HCs and nine SS patients (Supplementary 
Figure S2). 

MiR-21 (Figure 3a) showed upregulation in 10 out of 21 
samples (48%) (with an FC ranging from 1.4 to 4) (P<0.03). 
Interestingly, it resulted overexpressed (with an FC ranging 
from 1 .4 to 3) in five of five patients (100%) associated with a 
poor prognosis (mSS 30, 32, 39, 43 and 45), and at levels 
comparable to healthy donors in six of seven patients (86%) 
associated with a better outcome. These findings seem to 
indicate miR-21 as a more strictly potential prognostic marker 
in SS. In addition, as miR-21 maps at chromosome 17q23, a 
region frequently found amplified in SS,^ we tried to correlate 
the miR-21 expression with gains on 17q23 (highlighted with 
black bars, in Figure 3a), but we did not observe any 
significant correlation. miR-21 4 resulted over-expressed (with 
an FC ranging from 2 to 2000) in 17 of 22 cases (73%) 
(P< 0.004), whereas miR-486 was upregulated (with an FC 
ranging 1.8 to 29) in 18 of 23 patients (83%) (P< 0.0038). 
These data suggest that both these miRs represent useful 
diagnostic markers for SS. Moreover, interestingly, miR-21 4 
resulted absent or at very low level (comparable to HCs) in five 
of six long survivors (SS26, 36, 40, 51 , 59 and 63) analyzed by 
qRT-PCR, suggesting its active role in SS pathogenesis. 

l\/liR-21, miR-21 4 and miR-486 promote cell survival in 
Hut78 cell line. Because miR-21, miR-21 4 and miR-486 
overexpression were identified in this study as predictive of 
SS disease and associated to an unfavorable outcome, we 
investigated the potential effect on apoptosis and cell survival 
on the HUT78 CTCL cell line by modulating their expression. 
Hut78 shows miR-21 expression levels comparable to HCs, 
whereas it does not express miR-21 4 and miR-486 
(Figure 3). For these reasons, we followed an miR-21 loss- 
of-function and an miR-21 4 and miR-486 gain-of-function 
approaches. For miR-21 analysis, we transfected cells with 
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Figure 3 qRT-PCR validation of miR-21 (a), miR-214 and (b) miR-486 (c). mlRNA level was calculated for each SS pz relatively to a reference pool constituted by three 
HDs. In a, Black bars represented patients with a gain of 17q23 region, whereas white bars represented patients with normal ch 17q23 region. In gray are showed all other 
samples that were not analyzed by SNPs or CGH arrays 
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Figure 4 miR-21 , miR-214 and miR-486 increase apoptotic resistance in Hut78 cell line, (a) Apoptosis measured by nucleosome formation of HUT78 cells transfected with 
anti-miR21 and scrambled oligo, used as control, (b and c) Apoptosis assessed by flow cytometry 48 h post-transfection. Anti-miR-21 transfected cells showed a significant 
increase of apoptotic Annexin V +/PI + cells compared with control, (d) Apoptosis measured by nucleosome formation of HUT78 cells transfected with pre-miR-214 and pre- 
miR-486 compared with control (e) MTT assay showing increased cell viability of Hut 78 cells transfected with pre-miR-214 and pre-miR-486, compared with control. *P< 0.05. 



anti-miR-21 LNA-knockdown probe and a control-scrambled 
oligonucleotide to assess the biological effects 48 h post- 
transfection (n = 6). Transfection efficiency was >90% as 
determined by FACS analysis, and the resulting miR-21 
modulation was assessed by qRT-PCR (Supplementary 
Figure S3). MiR-21 knockdown significantly enhanced the 
apoptosis in transfected cells compared with controls as 
determined by the quantification of cytoplasmic histone- 
associated DNA fragments (Figure 4a) and by Annexin V/PI 
staining (Figures 4b and c). Conversely we did not observe a 
significant effect on cell proliferation as determined by MTT 
assay (not shown). 

For miR-21 4 and miR-486 studies, we transfected cells with 
pre-miR-214 and pre-miR-486 oligos and at 48 h post- 
transfection, we measured apoptosis (Figure 4D) and cell 
viability by MTT assay (Figure 4e) (n = 3). Results obtained 
revealed a decreased cell death associated to an increase 
of cell viability in miR-214 and miR-486 transfected cells 
compared with control. All these findings indicate that miR-21 , 
miR-214 and miR-486 contribute to the apoptotic resistance of 
CTCL cell line. 

Discussion 

miRNA whole gene analysis has proven to be a valid tool to 
identify diagnostic and disease progression signatures in 
several tumors and especially in hematological malignancies. 
In the present study, we have profiled for miRNA expression 
28 SS samples and six HCs, using microarray technology. A 
total of 45 miRNAs were the most significant among those 
differentially expressed between SS patients and HCs. 
Furthermore, we were able to identify 19 miRNAs with a 
potential ability of disease prediction. Until now, only one 
study has been done on miRNA expression for SS.^^ In this 
investigation, a set of 114 deregulated miRNAs able to 
distinguish SS from HCs was identified and 20 were the most 



discriminatory. Our results confirm, in part, these data 
showing a concordant increased expression of miR-214, 
miR-199a* and miR-7, along with the decrease of miR-342, 
miR-223, miR-92, miR-181a and miR-191. Conversely, a 
discordant expression level was observed for miR-145 and 
miR-1 8a. However, we have also found a large fraction of non- 
overlapping deregulated miRNAs that might be explained 
with the different methodological approaches used, such as 
miRNA platforms, the bioinformatics and statistical methods 
employed to analyze data and particularly to identify 
diagnostic and prognostic signatures. 

Our unsupervised analysis identified four sub-groups of SS 
patients and, thus, we wondered whether these different 
clusters were associated to any clinical parameters such as 
tumor progression, tumor burden, CD4/CD8 ratio, dominant 
TCR-V beta rearrangement, genetic lesions and survival. 
Interestingly, we observed a correlation with clinical outcome, 
recognizing a group of patients (belonging to cluster n.2) with 
the shortest survival respect to all other individuals enrolled in 
this study. We thus compared miRNA expression data with 
survival time using class comparison and risk prediction 
analysis. Following this approach, we obtained two indepen- 
dent lists of miRNAs differentially expressed between patients 
with unfavorable and favorable outcome (Table 1 and 2a). 
In the first one, miR-214 together with miR-1 99a and miR- 
199a* appeared the most significantly deregulated. The 
intersection of these two lists, identified 1 4 miRNAs (Figure 2). 
Among them, interestingly, we observed again miR-21 and 
miR-486 that we have studied in detail, miR-1 55 involved in a 
variety of tumors^° and miR-1 6 already found deregulated in 
SS."*^ Interestingly, many of these potentially prognostic 
miRNAs map to gain/loss genomic regions frequently found 
unbalanced in SS,^ suggesting that also genomic alterations 
might contribute to miRNA deregulation in SS. 

Several transcriptional or epigenetic mechanisms may 
influence miRNA changes. Narrowing our attention to the 19 
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predictive miRNAs, we noted that miRNA 21 4/1 99a cluster 
and 1 99a* are co- expressed in our sample set (Figure 1 b) and 
they were previously found to be upregulated by Twistl.^^ 
This is a transcription factor found overexpressed in SS,^ so 
that it might control the upregulation of miR-21 4/1 99a/1 99a* in 
this lymphoma. As a large fraction of our SS patients were 
already profiled for gene expression,^ we matched miRNA 
and mRNA results; in this manner, we noted that expression 
levels of miR-30c correlated with the expression of its 
host genes NFYC, suggesting a control at level of primary 
transcription rates. Similarly, miR-342 correlated with its host 
gene EVL that in turn is downregulated by 1 99a*. Moreover, 
we observed a number of tumor-related miRNAs like miR- 
18a, belonging to cluster 17-92, and miR-106b (cluster 106b- 
25) strongly involved in lymphoproliferative diseases. ^'^^ 
Additional members of the same clusters, such as miR-17- 
3p and miR-18a (cluster 17-92) and miR-93 (cluster 106b- 
25), are included in the list of prognostic miRNAs (Table 1 and 
2a). The genetic loci encoding the miR-17-92 (13q31) and 
106b-25 (7q22) clusters do not appear among the chromo- 
somal region recurrently altered in SS;^ therefore, their 
deregulation should be attributable to other mechanisms. 
For example, by c-Myc, found overexpressed in SS,^ that is 
able to induce the expression of these miRNAs in several 
tumors and animal models. 

In our dataset, we observed a diminished expression of let- 
7a, let-7b and let-7c (Figure lb and c). Let-7 family members 
are widely viewed as tumor suppressor miRNAs and they are 
strongly downregulated in many cancers and during tumor 
progression.^"^ Conversely, an upregulation of let-7d, let-7f 
and let-7i was mainly found in patients with unfavorable 
prognosis (Table 1 and 2a). A downregulation of miR-31 was 
observed in 22 out of 28 (79%) SS patients. miR-31 is an 
antimetastatic miRNA and its downregulation or deletion of its 
genomic locus at 9p21 (recurrently observed in SS) promotes 
invasion/metastasis of cancer.^^ These findings suggest 
that, similarly to other cancers, downregulation of miR-31 
might contribute to the atypical trafficking ability of these 
malignant lymphocytes, and might represent another potential 
diagnostic in SS. 

In this study, we focused our attention on miR-21 , miR-21 4 
and miR-486, because they most significantly correlated with 
diagnosis and prognosis of SS. Results indicated that miR- 
21 4 and miR-486 were over-expressed in the majority of SS 
patients therefore more likely related to diagnosis, whereas 
miR-21 upregulation is detected in a smaller number of 
patients thereby more probably associated to prognosis of 
SS. Moreover, performing functional experiments, we 
observed that all these miRs are able to promote cell survival 
in Hut78 cell line. Interestingly, miR-21 was also recently 
described to increase survival in another CTCL cell line,^^ 
whereas miR-486 and miR-21 4 have been recently asso- 
ciated to clinical outcome of many tumors. ^^'^^'^^'^^ Moreover, 
miR-21 4 over-expression has also been recently described 
to promote proliferation of healthy T lymphocytes via PTEN 
targeting.^^ Among the thousands of miR-21 gene targets, 
one of the best characterized and involved in neoplastic 
transformation is PTEN. Thus, it is possible that a concurrent 
over-expression of miR-21 and miR-21 4 increases SS cell 
survival, repressing the PTEN function. In conclusion, our 



findings indicate that miR-21, miR-21 4 and miR-486 might 
represent novel diagnostic/prognostic biomarkers for this 
incurable variant of CTCL. 

Materials and Methods 

Patients and healthy controls. Peripheral blood samples from 22 SS 
cases, six follow-ups of patients SS30, 33, 38, 39, 43 and six HCs were analyzed. All 
patients were enrolled in this study approved by the Ethical Committee of the Istituto 
Dermopatico dell'lmmacolata and informed consent was obtained in accordance 
with the Declaration of Helsinki. Diagnosis of SS was based on described criteria.^° 
The major clinical and immunological characteristics have been reported in 
Supplementary Table SI and obtained as described.^ 

miRNA expression profiling. Neoplastic and healthy lymphocytes were 
purified by dynabeads according to manufacturer's instructions (Invitrogen, Oslo, 
Norway) and total RNA extracted and quantified as already described in details.^ 
Labeled RNA was then hybridized on Agilent Human V1 miRNA microarray (Agilent 
Technologies, Santa Clara, CA, USA), consisting of 470 miRNA probes, according 
to the manufacturer's procedure. Agilent scanner and the Feature Extraction 10.5 
software (Agilent Technologies) were used to obtain the microarray raw-data. 

Bionformatics analysis. Data analysis was performed, by using customized 
R language-based script (www.r-project.org), using the Bioconductor packages 
(www.bioconductor.org). Quality control analysis, data normalization and 
hierarchical cluster were performed as described.^^ Data matrix of normalized 
values was subjected to class comparison and prediction analysis using BRB- 
ArrayTools version 3.9.0.^^ Class comparison allowed identifying the miRNAs 
differentially expressed between SS patients and HCs using a univariate 
two-sample f-test with a significance threshold of 0.01 , while controlling either the 
number or proportion of false discoveries. Class prediction was assessed to classify 
and to predict in which class was enclosed a sample in relation to its expression 
profile, using multivariate classification methods showed in Supplementary Table 
S4. Hierarchical clustering of the data was performed with the bioconductor 
heatmap package, using centered correlation and complete linkage similarity metric. 
Cluster Stab bioconductor package was used to compute cluster stability scores. In 
addition, we performed a survival risk prediction using BRB-ArrayTools version 
3.9.0.^^ This tool was used to directly predict the survival risk group based on 
expression data and two covariates as sex and age. MiR-21 target genes have been 
predicted applying mlR and an algorithm-associated MirBase software (http:// 
microrna.sanger.ac.uk). 

Survival and statistical analyses. A time-to-event analysis was performed 
using nonparametric KM product-limit survival estimates, and differences between 
KM survival curves were analyzed using the Wilcoxon test. Statistical analyses were 
performed using Statistical Package SPSS v.13.0. Results of qRT-PCR were 
expressed as mean plus or minus S.D. Statistical analysis was performed using the 
two-tailed Student f-test or one-way ANOVA. 

Quantitative real-time PGR (q-RT-PCR). qRT-PCR was used to confirm 
our microarray data. Selected mature miRNAs were analyzed using TaqMan 
microRNA assay (Applied Biosystems, Foster City, CA, USA). Z30 was used as an 
internal normalized reference. Quantitative analysis was performed by the ACT 
method as described.^^ 

Cell line transfections, cell viability and apoptosis. Hut-78 cell line 
established from peripheral blood of CTCL patient from American Type Culture Collection 
(TIB161) was grown in RPMI (Sigma Aldrich Biochemicals, St. Louis, MO, USA), 
supplemented with 10% heat-inactivated FCS and antibiotics at 37°C in a humidified air 
atmosphere at 5% CO2. LNA-knockdown miR-21 or LNA-Knockdown probe Scramble 
molecules 5'fluorescein (Exiqon, Vedbaek, Denmark), pre-miR-214 and miR486 or FAM 
Pre-miR-scrambled negative control (Ambion, Foster City, CA, USA) were transfected using 
the Amaxa Biosystems nucleofector electroporation system as described.^ Briefly, 2x10^ 
of Hut78 were resuspended in 100|il of nucleofection R, with oligos used at 30 nM. After 
electroporation, the cells were seeded into 24-well plates and the medium was refreshed 
~8h after transfection. The cell viabilrty and apoptosis were examined 48 h post- 
transfection with 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) reduction 
(Sigma -Aldrich Biochemicals) and quantifying cytoplasmic histone-associated DNA 
fragments (Roche Molecular Biochemicals, Manneheim, Germany) assays as described.^ 
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Flow cytometry. Flow cytometric analysis for apoptotic cells was performed by 
staining cells with FITC-Annexin V and propidium iodide, using the apoptosis 
detection kit (Clontech, Mountainview, CA, USA) Samples were acquired on a flow 
cytometer (FacsCalibur, BD Biosciences, San Jose, CA, USA), and analyzed using 
CellQuest software (BD Biosciences). 
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